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ABSTRACT 



Any discontinuity in flow at a plane boundary is dispersed by the diffusion of vorticity 
into a boundary layer and wake Experimental measurements in the past have been 
correlated with empirical formulations A generalization gives both the horizontal 
component and the vertical component of the mean velocity In the laminar sublayer 
the velocity is a solution of the diffusion equation In the turbulent boundary layer 
the velocity can be expressed by a Fourier integral In the free stream there is a 
vertical persistence of velocity The computation of velocity is provided by subroutines 
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INTRODUCTION 



The drag on a ship is partly a viscous resistance and partly a wave resistance The 
viscosity of the fluid causes a transport of momentum from the ship to the fluid The 
momentum is located in a boundary layer at the surface of the ship, and in the wake 
behind the ship The boundary layer is laminar at the bow, but the boundary layer is 
turbulent at the stern For a nonslip boundary condition the velocity relative to the 
ship is zero at the surface and is the velocity for free flow far from the surface A 
knowledge of the velocity may be found from an analysis of the analogous problem 
in the flow over a plane 

When a thin plate moves edgewise through a viscous fluid the plate entrains fluid 
if there is a nonslip boundary condition at the surface of the plate Each differential 
element of surface of the plate sets up a current which trails downstream and spreads 
out by diffusion The velocity at the plate is the accumulated sum of velocities in 
currents which have been created upstream 

The boundary layer of a plate has been the subject of many investigations 1-7 There 
are three principal methods of analysis The first method is statistical It leads to a 
system in which the number of variables is greater than the number of equations. 
Various schemes for closure have been proposed The second method is polynomial 
The velocity is expressed as a power polynomial in a limited region of a laminar 
boundary layer The third method is spectral The velocity is expressed as a Fourier 
integral The rate of change of the Fourier amplitude is determined by an 
mtegro-diflerential equation The evolution of a velocity distribution from an initial 
distribution is determined uniquely by the mtegro-diflerential equation 

The Fourier analysis has been applied to homogeneous isotropic turbulence Contact 
of the fluid with a solid boundary is the origin of real turbulence 

A comprehensive source of information about boundary layers is Boundary- Layer 
Theory by Schlichting 7 The Fourier analysis of flow in a fluid has been investigated 
by Batchelor 8 and by Orszag 9 



EQUATION OF CONTINUITY 

A fluid consists of particles with random velocities. The number of particles per 
unit volume determines a median density p Across any mathematical boundary there 
is a flux of particles from both sides A boundary which is moving at such a speed 
that the flux is zero defines a median velocity v 

In a continuous fluid the law of conservation of mass requires that within any 
mathematical boundary the mass density p and the particle velocity v are related by 
the equation 



dp w 

— dr 

at 



pv-ds = 0 



:d 



where t is the time, dr is a volume element within the boundary, and ds is a surface 
element on the boundary Application of the Gauss theorem leads to the equation 



dp 

at 



V-(pv) 



dp 

dt 



+ v-Vp pV-v = 



dp 

dt 



pV-v = 0 



( 2 ) 



The divergence of velocity is the rate of expansion of the fluid 
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EQUATION OF MOTION 

The bulk modulus k is the rate of decrease of pressure per unit rate of expansion 
A pressure pulse is propagated through the fluid with the speed of sound (tc p ) 1 
In a continuous fluid the gradient of velocity satisfies the identity 

Vv = |(Vv - V'v) 4 - |(Vv - V'v) (3) 

where V'v is the transpose of Vv The antisymmetric part is the rate of rotation Q and 
the symmetric part is the rate of strain 0 There would be spinup of fluid particles 
if stress were not symmetric In an isotropic fluid the pressure is isotropic and shear 
stress is proportional to shear rate The strain rate is partly an isotropic expansion 
rate, which is the trace of the strain rate, and is partly an anisotropic distortion rate, 
which is the shear rate The shear stress arises from the anisotropic distortion rate 
The stress Z is given by the equation 

Z = - p I ~ §M I + 2/j, 0 (4) 

where p is the pressure, and p is the viscosity. 

The force on any mathematical boundary is given by the Gauss theorem 

j V-E dr = J E ds (5) 

where dr is a volume element within the boundary and ds is a surface element at the 
boundary The equation of motion per unit volume is given by the equation 

dv 

P^=v-E (6) 



As a consequence of the identity 



V-V'v = VV-v 



(?) 



the equation of motion is the Navier-Stokes equation 



dv 

— 4- v- Vv = - V 

dt 



— + 5 - VV-v + - V-Vv 
P P P 



The kinematic viscosity v is given by the equation 



( 8 ) 



v = — 

P 

For hydrodynamics the kinematic viscosity is constant, 
equation of motion leads to the equation 



dt 



(V-v) ■+• v-V(V-v) 4 - (Vv-Vv) = - V-V 



if 



( 9 ) 

Differentiation throughout the 

+ i - V-V(V-v) (10) 

P 



where a tensor enclosed in parentheses is the contraction of the tensor by internal 
scalar multiplication. This equation is a Poisson equation for the determination of 
pressure Differentiation throughout the equation of motion leads to the equation 



— (Vxv) 4 * v-V(Vxv) + [Vv-Vv] = — V-V(Vxv) 
dt p 



(ID 



where a tensor enclosed in brackets is the contraction of the tensor by internal vector 
multiplication This equation is a diffusion equation for the determination of vorticity 
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BLASIUS PROFILE 



For two-dimensional flow it is fashionable to define Cartesian coordinates x, y, 2 
such that x is parallel to the plane in the direction of flow, y is perpendicular to the 
plane, and 2 is parallel to the plane in the direction perpendicular to the flow Let 
i,j,k be unit vectors in the directions of increasing x, y, 2 The particle velocity v is 
given by the equation 



v = ui + v j (12) 

where the components u, v are zero at the plane but approach U, 0 with increasing y 
Continuity in an incompressible fluid is expressed by the equation 



du dv 
dx dy 



(13) 



The continuity equation is satisfied when the components of velocity are derived from 
a stream function y in accordance with the equations 



u = 



dy 

dy 




(14) 



This follows from the principle that a second-order derivative is independent of the 
order of differentiation 

For stationary flow the derivative with respect to t is zero Far downstream derivatives 
with respect to x approach zero, and the Navier-Stokes equation is reduced to the 
equation 



du du d z u 

4- v - V T 

dx dy dy z 



15) 



Far downstream the profile of the boundary layer approaches a constant limit 

The limiting profile is given by a relation between two dimensionless variables 
f and <f, which are defined by the equations 



£ = 




y 

Vx 



y = V vUx ( (16) 



Differentiation leads to the equations 



u = U(' 



v - 



vU , 

— (W - 0 

x 



(17) 



Then substitution in the reduced equation of motion leads to the Blasius equation 

<rr + 2<r = o (is) 

Boundary conditions for this differential equation are given by the equations 

<■ = <" = 0 at £ = 0 (19) 

- 1 at £ - 00 (20) 

In the original derivation by Blasius 1 the variables were so defined that (' approached 
a limit of 2 The derivation herein follows the derivation by Schlichting 7 where the 
variables are so defined that (' approaches a limit of unity. 
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An ascending series is given by the equation 



C= I (2D 

m = 0 

This series satisfies the boundary conditions at £ = 0 Substitution in the differential 
equation leads to the recurrence equation 

m- 1 

£ (3m - 3k - 2)(3m - 3k - l)a t a m _ t _, + 2(3m)(3m + l)(3m + 2)a m = 0 (22) 

fc = 0 

The recurrence is started with a value of a 0 which is adjusted by trial to make the 
series meet the boundary conditions at £-*<». The computed value of a 0 is 0 16603 
The ratio between successive coefficients of the series tends to a constant limit with 
increasing order It is possible to estimate a remainder after a finite number of terms 
on the basis of the geometric series. In no case can the series be used for £ greater 
than 5.69, where the ratio between terms becomes unity 

For large values of £ the value of is the sum of a linear term and a correction. 
Correct to first order the correction is a solution of the equation 



(£^cR" + 2C'" = 0 (23) 

where c is an arbitrary constant. The computed value of c is - 1 72077 The solution 
of the differential equation is given by the equations 

r* OO 

c = $ + c -4(£ + c) J e~? (?+c)2 d* + ZA e'i (c+c)2 (24) 

r% oo 

C'=l-^J e'* u + c)Z d$ (25) 

where A is an arbitrary constant The arbitrary constant is selected so as to make 
the corrected linear terms coincide with the series expansion at £ = 5 0 The integral 
in the correction is given by the equation 

f »oo 

J e'i (?+e)2 a!£ = Vn [l - erf(±(£ + c))] (26) 

where erf is the error function. 

Previous analyses and computations on the Blasius profile are summarized by the 
data in a table in a paper by Howarth 2 , which is the basis for the Table 7 1 in the 
text by Schlichting 7 . That table of data is reproduced by the following subroutine 

SUBROUTINE BLSSPF (AU. AN, AX, AY, FU, FV) 



FORTRAN SUBROUTINE FOR BLASIUS PROFILE 
********************************************************************************************* 

The free-stream velocity U is given in argument AU, and the kinematic viscosity v 
is given in argument AN. The coordinates x, y are given in the arguments AX, AY Series 
expansions and error integrations are used in the evaluation of the Blasius profile. 
The components u, v are stored in functions FU, FV. 
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TURBULENT PROFILE 



Far downstream the flow in the boundary layer is turbulent There have been many 
experiments on the turbulent boundary layer In the absence of a fundamental theory 
the experiments have been correlated with empirical formulations A celebrated 
formulation is the logarithmic law, which expresses mean velocity as the logarithm 
of the distance from a wall The logarithm cannot be used at the wall where velocity 
is zero and the logarithm is -«>, or at infinite distance from the wall where velocity 
is finite and the logarithm is +«>. Corrections have been published by Reichardt 4 and 
by Thompson 5 The argument of the logarithm is incremented by unity and the logarithm 
is blended with the free-stream velocity 

The Reynolds number R x is defined by the equation 

Ux 

R x = — (27) 

v 



where U is the free-stream velocity, v is the kinematic viscosity and x is the distance 
downstream. The drag D on the wall is given by the equation 



D = ± pU z 



cx 

(log*,) 2 - 58 



(28) 



where the constant c is defined by the equation 

c = 0.455 (log 10) 2 58 (29) 

This formulation is equivalent to Equation 21 16 in the text by Schlichting 7 A shear 
velocity u* is defined by the equation 

u* = J- (30) 

\p 



where r is the shear stress at the wall, and p is the density The shear stress r is 
given by the equation 



dD_ 

dx 



(31) 



The boundary layer thickness 6 is determined by the equation 



u+6 _ R x 
— = (log*,) 2 - 58 



(32) 



This equation is equivalent to Table 21 1 in the text by Schlichting 7 Finally the mean 
velocity u is given by the equation 



yxu 

u = log 

tC 




+ (1 -7)i/ 



(33) 



where /c is a constant and y is the blending function The experimental value of tc is 
given by the equation 

tc = 0.40 (34) 



Various empirical schemes for the blending function have been tried, but a more 
logical basis for the blending function would be the error function To within 
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experimental error the blending function can be expressed by the equations 



y 

n = n 6 



7 = 



, e rf(r) + f ) - erf(7? - f ) 
erf(f) 



( 35 ) 



Then the blending function is in the range 1 > y > 0 while y is in the range 0 < y < °° 
Differentiation with respect to x leads to the equation 



_ i 



pU z 



(log R x ) z 58 

and further differentiation leads to the equation 

2.58c 



1 - 



2.58 
log R x 



dr 

dx 



= - ip u 3 



x(log R x ) 2 



1 - 



3.58 

log/?. 



Then the derivative of u* is given by the equation 



du* 



1 dr 



dx 2 * t dx 



Differentiation with respect to x leads to the equation 



U 



dt 5 



dx u*(\ogR x ) 



3.58 



1 - 



log R x - i(2.58)- 



1 58 
log R x 



1 - 



2.58 
log R x 



Then the derivative of y is given by the equation 



— = V n - - 

dx 6 



/ TT v 2 / TT . 2 

‘(*1 - jr ) -( 17 +- 5 -) 



erf(|) 



2 l d6_ 
6 dx 



(36) 



(37) 



(38) 



(39) 



(40) 



Finally the derivative of u with respect to x is given by the equation 



du 

dx 




1 + tc 



u*y 



d.y 

u. — + y 
dx 



du * 
dx 



u*y 

v 



1 + /c 



du. dy 

7— U — 

dx 



u*y dx 



(41) 



Although the derivatives are given by empirical functions with finite numbers of terms 
their integration with respect to y would lead to an unlimited number of terms. In a 
practical integration with respect to y the integrand is best approximated by a power 
polynomial in y and integrated term by term. A discrete set of integrands is converted 
into a discrete set of coefficients by 11 -point Lagrange interpolation and the integration 
is completed coefficient by coefficient. 

The stream function if/ is given by the equation 

V' = I u dy (42) 

Jo 
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and the vertical component v is given by the equation 

dy 



P ^ 
Jo dx 



( 43 ) 



The integrations are performed with the aid of the following subroutine 

SUBROUTINE T8LNP r (AU, AN, AX, AY, Fu, FV) 



********* 






********************* 



FORTRAN SUBROUTINE FOR TURBULENT PROFILE 



k****************** 



********* 



k*************** 



********************4 



The free-stream velocity U is given in argument AU, and the kinematic viscosity v 
is given in argument AN The coordinates x,y are given in the arguments AX.Ar 
Empirical formulations and Lagrange interpolation are used in the evaluation of the 
turbulent profile The components u, v are stored in the functions FU, FV. 

LINE FLUX PROFILE 



The usual assumption is made that pressure is constant over a fiat plate in a steady 
flow parallel to the plate Let the free-stream velocity + £/i be disturbed by a small 
counter velocity v. The Fourier transform can be applied to a function which is zero 
everywhere except at the origin Then Fourier integration leads to the equation 

v= -^LUU) c '''*“ vd ' c (44> 



That this expression for velocity has zero divergence can be verified directly by 
differentiation Furthermore it can be derived from the stream function ip in the 
equation 



V = ~ 



2?t J.oo ifc 



(45) 



At x = 0 an application of the Euler theorem and an integration of the sine quotient 
function shows that ip = - | for all y > 0 and ^ | for all y < 0. Thus the flux in the 

counter flow is unity at x - 0 The curl of the velocity is given by the equation 

I /» + oo / \ 2 2\ 

Vxv =4- — — — e - x * + **v dx = yk (46) 

2tt J_ 00 itc 



To within small quantities of second order the vorticity satisfies the differential equation 



dy 

U— = v 
dx 



d^y 

dx 2 




(47) 



Substitution of the Fourier integral into the differential equation leads to the equation 

- UX = v (X 2 - k 2 ) (48) 



which may be solved by the quadratic rule to give the equation 

^ _ _ U_ ^ \( U \ 2 2 2vk 2 

X ~ 2u + ^\2u) +K " u + y/l T z + 4 



(49) 



The parameter A is an even function of the parameter fc The radical is negative 
upstream and is positive downstream Otherwise the integrals do not converge 
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When y is increased to infinity the integration with respect to tc makes a significant 
contribution to the integration only where k - 0 
Upstream the parameter A is given by the limit 

A - - — (« - 0) (50) 

1 / 



Thus the stream function ^ 1S given by the limit 




and the velocity v is given by the limit 




(y - ± 00 ) (5i) 



(y- ±~) (52) 



The vertical component of velocity persists with increasing y, but diminishes rapidly 
with distance |x| upstream. 

Downstream the parameter A is given by the limit 

A - 0 (tc - 0) (53) 

Thus the stream function ip is given by the limit 

- * 5 (y - ± “) (54) 

and the flux in the counter current is everywhere unity. It is independent of distance 
x downstream. 

In the integrals for the components of velocity the integrands are the products of 
the monotonic factors 




A 

2m*; 



e 



~\x 



(55) 



and the oscillatory factor 



e~ y (56) 

The integration through any number of cycles of the oscillatory factor can be completed 
if the monotonic factors are expressed as power series in tc. The range of the power 
series is limited by the presence of the radical 



U \ Z * 
— + /c 2 

2 v) 



(57) 



in the parameter A For small values of tc the radical is expressed by the equation 

v 2m 

(58) 



U 



Zv 



8 = U_ y (- l) m *‘(2m)! !2wc\ z 
2v^{2m- l)2 am (m!) 2 l U ) 



This series converges only when tc meets the limitation 

U U 

< tc < + — 

2v 2v 



(59) 



The series is an even function of tc. For large values of tc the radical is expressed by 
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the equation 



'/_tn 2 + 2 = “ (-ir‘(2m)' / u \ 2m 

\2v) +K K m t 0 (2m-l)2 2m (m'H2w/ 



This series converges only when tc meets the limitation 



2v 1 2v 

— < - < + — 

U tC U 



(60) 



(61) 



The series is an odd function of 1 / tc Efficient evaluation of the series is only possible 
if the ranges of their arguments are much less than the limits of their convergence 
The monotonic factors are expanded in each of a sequence of intervals of limited 
range 

The value of tc in the first interval is given by the equation 

tc = t) + 6 (62) 

where 77 is the center of expansion and 9 is the variable of expansion The variable 9 
is given by the equation 

9 = eu (-1 £ u £ +1) (63) 



where e is half the range of expansion and u is a variable of interpolation The first 
interval straddles the origin where 77 = 0 , then in subsequent intervals 77 is incremented 
by 2e Thus the monotonic factors are approximated by the series 

CO ^ 00 

= I a m u m -e’ Xx = £ a m u m (64) 

m = 0 m = 0 



Required for the computation of velocity are the integrals 

/ m du 



e e 



ir?y 



r 

Jo 






e e 



XT)y 



a m u m e xcyu du 



Required for the integration is the recurrence equation 

*±i 

u m e xcyu du = 

f 0 



f 

Jo 



u 



ley 



icyu 



±l 771 f ±1 

u m_1 e xcyu du 

o tey Jo 



The recurrence is started with the initial integral in the equation 



f 

J o 



du = 



e xcyu 7tl 

Jo 



The recurrence is cycled in ascending order if ey satisfies the inequality 

|ey| g 17 



(65) 



( 66 ) 



(67) 



( 68 ) 



Otherwise the recurrence is cycled in descending order with an initial approximation 
for 77 i = 64 

The value of tc in the last interval is given by the equation 



tc = 




(-lSu^tl) (69) 



where c5 is the limit of integration The last interval straddles the point where £ = 0 
Continuity through the point is achieved by giving to the radical in X the same sign 
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as the sign of < The monotonic factors are approximated by the series 






-(X - K )x _ 



= E C m u r 



-e- (x - K)x = E c m u r 



m - 0 m = 0 

Then the components of velocity are expressed in terms of the integrals 



oo r* a 

V C 

L m 

m =0 v 6 



- e ~ K,x ~ t v > die 



Required for the integration is the recurrence equation 



_L 2 «(x iy) ^ - - 

/c m 



e - k(x-i y) 

(m - l)/c m ‘ 



(j ~ iy) r 1 

(m - 1) 



f — 1 — e 

J a 



I -' y) die 



The recurrence is started with the initial integrals in the equations 

-<(x-xy) -|°° 

€>-<(x-vy) dK = 



* - iy h 



J i dK = J 



-< 5 (x-xy) « 



— dt = - Ei(-6(x - iy)) 



(70) 

(71) 

(72) 

(73) 

(74) 



The recurrence is cycled in descending order if 6{x - iy) satisfies the inequality 

|<5(x - iy)| S 17 (75) 



Otherwise the recurrence is cycled in ascending order with an initial approximation 
for m — 64. 

For the series expansions the arguments have Chebyshev spacing and interpolations 
are made with 11 -point central Lagrange interpolation. Preliminary computations have 
established a matrix of coefficients such that the coefficients for progressively increasing 
powers of the argument are obtained with the product of an array of values of the 
function and the matrix of coefficients. 

If x satisfies the inequality 



7 u [nay u la 

8u > \\ x ) + v x 



(76) 



then e and 6 are given by the equations 

6 



6 = 




U IQ 



v x 



(77) 



and the last interval can be jettisoned without significant error Otherwise € and 6 
are given by the equations 



U_ 

Qv 




(78) 



in which case the Reynolds number is less than 36. 

The components of velocity for the line flux are computed by the following subroutine. 
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SUBROUTINE LNFXPF (AU, AN, AX. AY, FU. FVJ 



FORTRAN SUBROUTINE FOR LINE FLUX PROFILE 
********************************************************************************************* 

The free-stream velocity U is given in argument AU, and the kinematic viscosity v is 
given in argument AN. The coordinates x, y are given in the arguments A/, Ax Ascending 
and descending recurrence relations are used in the evaluation of Fourier integrals 
for a line flux profile The components u,v of velocity are stored in functions Fu. Fv 

NONSLIP BOUNDARY 

Exploratory computations with LNFXPF have indicated the nature of the flow from 
a line flux 

If x and y are decreased to zero, the contribution to integration extends to large 
values of k where A approaches the approximation 

* - - — * kl (kl - °°) (T9) 

At the limit of small x and y the velocity is given by the equation 

1 ! x|i + yj 

v = 2 (x - 0, y - 0) (80) 

TT x d + y d 

which is an efflux upstream and an influx downstream 

When x is increased to infinity downstream, the integration with respect to k makes 
a significant contribution to the integration only where /c - 0 and A approaches the 
approximation 



A - 




(fc - 0) (81) 



Then the velocity is given by the equation 



v = 



U 

47T vx 




£ 4ux 



(x-oo) (82) 



This is a Gaussian distribution of velocity 

In a continuous distribution of line fluxes the efflux and the influx at the leading 
edge would tend to cancel It is assumed that for integrated distributions the Gaussian 
profile is adequate 

The formula for velocity does not meet a nonslip boundary condition because the 
velocity at y - 0 varies inversely as the square root of distance downstream However, 
a continuous distribution of line fluxes may be integrated to give a constant velocity 
As a consequence of the equation 



f 

Jo 



da 
y/a vx 



= 2 



a 

\U - a 



= 7T 



JO 



(83) 



the integrated velocity is given by the equation 




da 

v(x - a) — — = - Ui 
v a 



(y = 0) (84) 
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as required by the nonslip boundary condition Thus the strength of line flux per unit 
displacement is given downstream by the expression 






4 vU 
rra 



(85) 



The Gaussian distribution is equivalent to the dispersion of a line pulse of vorticity 
which has diffused outward for a time equal to the expression 



x - a 
U 



(86) 



and has been swept downstream with the speed U 

FOURIER TRANSFORM 

Let f(x) be a periodic function of x. Then it can be approximated by the Fourier 
series 



l a o + E a n cos ni+ J 6 n sm nx 



(87) 



Inasmuch as cosines are even and sines are odd, they satisfy the equation 



r 



cos kx sin mx cLx = 0 
Application of the Euler theorem to the following integrals 

x(k+m)x dx ± 



(80) 



J + rr p + n 

e i(*+m)x dx ± i e t(fc ’ m)x 



dx - | 



" e x(*+m)x - 


+• TT 

+ i 


' g x(k-m)x - 


_ i(k + m) _ 


± 2 

- rr 


1 

? 

1 

-V 



= 0 (k * m) (89) 



shows that the trigonometric functions satisfy the orthogonality relations 



p + f! 

J -n 



cos kx cos mx dx 



-r 

0 - n 



sin kx sin mx dx = 0 



Application of the addition theorem leads to the equations 



I 



dx - 2rr 



p + rr p + n 

cos 2 kx dx = 

J - rr J - n 



sin 2 mx dx = n 



The mean square error for the Fourier series is given by the equation 

J r> + 7T OO OO 

a 2 = — )^a 0 + £ a n cos m + £ 6 n sin nx - f(x) \ z dx 

' ^ J -n n= 1 rx= 1 



(90) 

(91) 

(92) 

(93) 



Differentiation with respect to the coefficients shows that for least squares error the 
coefficients are given by the equations 



i r w 

a n = - 

77 J-n 



f{x) cos nx dx 



- 1 r n 

77 J-n 



f(x) sin nx dx 



(94) 



Because of the symmetry of the trigonometric functions, the trigonometric series is 
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given by the equation 

n = + oo /» + tt 

/(* 

The substitutions 



. n= +oo /» + 

-s.lL 



. n = +«> p + rr 



/(s) COS n(x - s) ds 



f(s) e 



tn ( x - s) 



ds 



2ttx 



2ns 



s - 



(95) 



(96) 



expand the range of approximation from 2n to L The functions f(x) and f(s) are 
replaced by the functions F(x) and F(s) as expressed by the equation 



F(x) 



= - U E f 5 

l n—J-i 



Zirn , . 

F(s) e L ds 



The substitution 



2nn 



k = 



(97) 



(98) 



replaces summation with respect to n by integration with respect to k in the limit as 
L - «. The approximation of summation by integration requires the functions to have 
properties of integrability The Fourier series is approximated by the Fourier transform 



F(x)e- 



r> + oo 

F(x) = A(fc) e lJCX dtc 

fj — oo 



dx 



(99) 



(100) 



In a multidimensional space there is a Fourier transform for each coordinate Each 
transform in a series of transforms can be applied to the amplitude of the previous 
transform 

The two-dimensional transform is given by the equations 



A(a, (3) = J J f( x , y) e l<a *" tfy) dx dy 



F(x, y) = A(a, (3) e liax * fiv) da df3 



( 101 ) 
( 102 ) 

Let <p(x,y,z) be defined by the equation 

cp(x, y, z) = j J /l(a,«e-^^ 2|t|t,(aitflv) dad/S (103) 

Then ^ is a solution of Laplace’s equation 

d z <x d z <x d z <p , 

v-v^ = — | + -f- + — ? = 0 (104) 

dx dy dz 

wherever z * 0 The derivative of <p with respect to 2 is given by the equation 

_ g = ± J J v'a* + /3 2 A(a. 0) da d p (105) 

where the sign is - for z < 0 and the sign is + for z > 0 The difference in the derivative 
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on opposite sides of the plane 2 = 0 is 4 ncr where a is the source density on the plane 
The amplitude 

2 Va 2 - (3 Z A(a. (3) (106) 

is the amplitude for 4 ncr For a unit source at the origin the amplitude is given by 
the equation 



4(a, p) = 



1 



2 tt va 2 



107) 



Two solutions of Laplace's equation are identical to within an additive constant if they 
have the same normal derivative on a boundary Thus the inverse of distance is given 
by the equation 



1 



V J 2 + y 2 



-i- 2 2 2rr J j Va 2 -1- 



rNj a 2 ^ 2 |*|n(ax^y) da ^ 



(108) 



Application of the Fourier transform to the real exponential factor in the integrand 
leads to the equation 



; = i J J 7 ^775 



(109) 



The evaluation of this Fourier integral can be completed after a transformation of 
coordinates. 

Let /c, 0. 0 be polar coordinates with polar axis in the direction of r Then a,{3,y 
are given by the equations 



a = te sin 0 cos 0 (2 = tc sin 9 sin <f> 

da dp dy -» /c 2 sin 9 dtc d9 d<p 

and integration with respect to 9, 0 leads to the equation 



y = tc cos 9 (110) 

( 111 ) 



sin ter 



dtc 



ter 



( 112 ) 



This equation confirms the existence of a Fourier transform for the potential of a 
pole 

Let r, 9 , 0 be polar coordinates with polar axis in the direction of te. Then x, y , 2 are 
given by the equations 

x = rsin0cos0 y = rsin0sin0 z = rcos0 (113) 

dx dy dz - r z sin 9 dr d9 d<p (114) 

Another solution of Laplace’s equation is given by the equation 

d I 1 \ cos 9 

dz \ r / r 2 

Its Fourier amplitude is given by the equation 



(115) 



A{tc) = J j j ) e lKTC09e r z sin 9 dr d9 d<p (116) 



Integration with respect to 9, 0 is completed with an integration by parts to give the 
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equation 



4(0 



i 

27T 2 AC 




cos(/cr) 



sin (ter) ) dr 
ter ) r 



(117) 



Integration with respect to r is completed with an integration by parts to give the 
equation 



4(0 



l 

2n z te 



sin(/cr) 

Jo 



i 

2tt 2 /c 



(118) 



This equation confirms the existence of a Fourier transform for the potential gradient 
of a pole. 

The gradient of inverse distance is given by the equation 



xi + yj + zk 
\x z + y 2 + z 2 ] 2 




a i + /Sj + y k 
a 2 + /3 2 + y z 



e l(ar ^ y+n ) da dp dy 



(119) 



The velocity of a point vortex of unit strength is given by the equation 



v = 



— kxV 
4tt 



r 



_L_ 

47T 



yi - *j 



3 

2 




ffi ~ Qj 

a 2 + £ 2 + y 2 



e x(«***V + 7«) dad 0 



( 120 ) 



This equation confirms the existence of a Fourier transform for the velocity of a point 
vortex. 

The circulation around a circle at coordinates r, 6 is given by the equation 



I 



vdr 




( 121 ) 



A change from spherical polar coordinates to cylindrical polar coordinates and 
integration with respect to z with the aid of the equation 

+ y 2 dz 
z z Vx 2 + y 2 +z 2 

confirms that the circulation is unity around a line vortex of unit strength and of 
infinite length. 




Vi 2 + y z 



,2-] + 



( 1 



INCOMPRESSIBLE FLUID 

Let r be a position vector in physical space and let k be a position vector in wave 
number space. The position vectors are defined by the equations 

r = xi + yj + zk k = ai + {3} 4- yk (123) 

where i,j,k are orthogonal unit vectors 

Let the velocity vector v at the position vector r be expressed by the equation 

v(r) = J A(k) e lt r |dK| (124) 

where k is a vector in wave number space, jd>c| is a volume element in wave number 
space, and A(k) is the amplitude of the Fourier component with wave number k The 
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amplitude is given by the equation 



^3 j v < r) *' 



the equation 



V- v = i J k-A(k) e xx r 



and the curl of velocity is given by the equation 

Vxv = i J kxA(k) e ix 
The Laplacian of velocity is given by the equation 

• 2a/v' 



r r |dr| 


(125) 


The gradient of velocity 


is given by 


d kI 
n 


(126) 


' \cLk\ 


(127) 


r |dK| 


(128) 


: ' r |dK| 


(129) 



The amplitude k-A(k) of the divergence V-v is given directly by the equation 



k-A(k) = — I V-ve r |dr| 
87T j J 



(130) 



If the divergence is zero for every r in physical space, the amplitude is zero for every k 
in wave number space. The amplitude satisfies the orthogonality equation 

k-A(k) = 0 (131) 



If a nonzero divergence did happen to occur in an incompressible fluid, a pressure 
pulse would be created, and the nonzero divergence would be dispersed by a potential 
flow 

Terms which are quadratic in velocity are expressed by double integrals Let K t and k 2 
be variables of integration in the double integration. Incompressibility is expressed 
by the equations 

k 1 A(k 1 ) = 0 (kj - k 2 )*A(k 1 - k 2 ) = 0 k 2 A(k 2 )=0 (132) 



When terms with the same wave number are collected in the integration the amplitudes 
are collected in a convolution. The product of velocity and its gradient is given by the 
equation 



v-Vv = i j J A(k 2 )k 1 -A(k 1 - k 2 ) e XX{ r IchcJIdK^ 
The square of the gradient is given by the equation 

(Vv-Vv) = - J J A(k 2 ) k 1 k,-A(k 1 - k 2 ) e " 1 r IctaJIdKa 
This scalar is the divergence of the vector 

i f f A(K 2 )-^ i ^ i -A(K 1 -K a )e" 1 ' r |dK 1 ||dK 2 | 

J J K r K i 



(133) 

(134) 



(135) 
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The gradient of pressure is given by the equation 



Vp = -pi 



// 



A(k 2 )- 



K, K, K 



K, K 



•A(kj - k 2 ) e l%l r IdKjlldKgi 



(136) 



Substitution in the equation of motion leads to the mtegro-diflerential equation 



which expresses the evolution of the spectral representation of the velocity in an 
unbounded fluid 

The rate of change of amplitude is expressed as the sum of three terms The first 
term replaces x in the Fourier transform with x - Ut The second term expresses the 
influence of amplitude at other points in wave number space The third term gives 
the rate of viscous dissipation 

The rate of change of A(kj) is influenced by the presence of A(k 2 ) in the vicinity 
of Kj In the integrand the postmultiplication of A(k 2 ) by a tensor eliminates any 
component of A(k 2 ) in the direction of Kj Other terms outside the integral are the 
products of scalars and A(kJ Thus the mtegro-diflerential equation preserves the 
orthogonality of Kj and A(kj) 

There is no contribution to the integration where k 2 is orthogonal to k 1 and A(k 2 ) 
is parallel to There is no contribution to the integration where k 2 is collinear 
with k j and A(Kj - k 2 ) is orthogonal to Kj If k 2 is on the perpendicular bisector of 
then interchange of k 2 and Kj - k 2 leaves amplitudes the same but reverses the sign 
There is no contribution to the integration by integration along the perpendicular 
bisector of Kj There is a maximum contribution to the integration when k 2 is on the 
perpendicular to Kj through the tip of Kj The influence of A(k 2 ) is a pattern which 
is crossed by nodal lines 



During unbounded evolution the velocity deviates gradually from the boundary 
conditions at the surface of a plate Velocity is injected gradually into the stream to 
maintain the boundary conditions It is only at the surface of the plate that velocity 
is injected Everywhere else the evolution of velocity is free Let the plane with plate 
be divided into an equally spaced grid Over each grid point there is a sine quotient 
function Each sine quotient function is unity at its own grid point and is zero on 
every grid line which does not pass through the grid point An analytic function can 
be expressed as a series in sine quotient functions The coefficients of the terms in 
the series are just the values of the function at the grid points The terms for defect 
in^ velocity are finite only for grid points within the area of the plate The sine quotient 
function at any grid point can be expressed by a Fourier transform The amplitude 
of the Fourier transform is a rectangle of constant density in wave number space 
Let x. y , 2 be Cartesian coordinates relative to a grid point with x positive downstream, 
with y positive across the stream, and with z positive perpendicular to the stream 
Let Q(a, (3) be the Fourier amplitude for the sine quotient function for the grid point 
and let A be the defect of velocity at the grid point The velocity is injected into the 
laminar sublayer where the absence of convection reduces the equation of motion to 



r; A(k 1 ) = - i U i-Kj A(kj ) - i 



dt 




DIFFUSION 
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the diffusion equation 



< 3 v / d z v d z v d 2 v\ 

dt ' dx z dy z dz z j 



( 138 ) 



The defect in boundary condition is concentrated initially against the plate surface 
and then diffuses outward The velocity for a defect A is given by the equation 




V4t /t 



Q(a. /?) da d p 



(139) 



in which the velocity is proportional to the complementary error function That this 
is a solution of the diffusion equation can be verified by direct substitution if z * 0 
The error function is symmetric with respect to the plane 2 = 0 The Fourier 
amplitude of the error function is the sum of an integration in the range -» < 2 < 0 
and an integration in the range 0 < 2 < +00 The integrals in the two ranges are complex 
conjugates Integration by parts in each range and cancellation of complex conjugates 
lead to the equation 

2 f 00 2 e~ yZvt 2 r+'y'tt z 

-^ 7 = e u dudz = 7 = e~ u du (140) 

Vrr J-I2 L i y Vtt J -xyVUi 




The amplitude of the error function is an even real function of 7 The velocity is given 
by the equation 



v = 



at r f 2 

2 n J J J vn i y 



e 



da d p dy 



(141) 



To verify this equation the error function can be replaced by its absolutely convergent 
power series, then term-by-term integration leads to the series for the arcsine of 
unity. The equation gives the initial velocity for the defect A The velocity continues 
to evolve thereafter in accordance with the full equation of motion for unbounded flow. 

DISCUSSION 



An analysis of the complete Navier-Stokes equation for Reynolds numbers less than 
1218 has been given by Schwiderski and Lugt 3 . They found overshoot in the velocity 
profile and the overshoot increased with Reynolds number. There does not seem to 
be any confirmation for the overshoot. There is no overshoot in the Blasius profile. 
Validation of the Blasius theory for Reynolds numbers more than 100000 is given by 
the excellent agreement between measurement and theory in Figure 7 9 in the text 
by Schlichting 7 It would be nice to have more experimental data for a smaller Reynolds 
number 

The flux of momentum across any cross section is equal to the drag upstream of 
the cross section The flux of momentum is given by the integral 




u(U - u) dy 



(142) 



where u is the instantaneous velocity in the fluid. The instantaneous velocity can be 
replaced by mean velocity only in the case of laminar flow. Otherwise the mean square 
of the instantaneous velocity is greater than the square of the mean velocity by the 
mean square of the fluctuation in velocity 
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Heretofore it has been assumed that the blending function for a turbulent boundary 
layer is symmetric with respect to the middle of the boundary layer, and is zero 
outside the thickness of the boundary layer It is more likely that the blending function 
diminishes gradually with distance An empirical formulation is compared with 
experimental data in the following figure 




6 

where the curve represents Equation (35) and the circles are from a tabulation in the 
report by Thompson 5 

The assumption that the horizontal component of velocity is constant everywhere 
in the free stream leads to a vertical component which persists to infinite distance 
from the boundary However, the distribution of vertical velocity initiates a jet and 
diffusion of vorticity disperses the jet in the free stream 

At the leading edge of a boundary layer the flow is laminar, but the flow becomes 
unstable where the Reynolds number is 520, and the flow becomes fully turbulent 
where the Reynolds number is 500000 For a free-stream flow at 20 knots the laminar 
flow is unstable at only one twentieth of a millimeter downstream, and the flow is 
fully turbulent at five centimeters downstream Any laminar flow in the boundary layer 
of a full-scale ship is insignificant 

CONCLUSION 

It is concluded that the most practical representation of the mean velocity is a 
generalization of the empirical formulation by Schlichting 
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